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).  Abstract 

A  review  is  undertaken  of  the  various  forms  that  have  been 
obtained  for  the  recurrence  relation  from  which  the  eigenvalues 
of  Laplace's  tidal  equation  may  be  obtained.  Such  forms  are 
shown  to  be  analytically  consistent  and  are  discussed  in  relation 
to  their  subsequent  numerical  evaluation.  By  determining 
eigenvalues  of  equivalent  depth  for  a  given  frequency  of 
oscillation  instead  of  the  other  way  around  the  problem  becomes 
a  straightforward  one  of  matrix  diagonalization.  If  solutions 
are  based  on  normalized  Legendre  polynomials  the  matrix  is 
symmetric.  A  method  of  evaluating  the  related  wind  functions 
is  described. 
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Abstract 

A  review  is  undertaken  of  the  various  forms  that  have 
been  obtained  for  the  recurrence  relation  from  which  the 
eigenvalues  of  Laplace's  tidal  equation  may  be  obtained. 
Such  forms  are  shown  to  be  analytically  consistent  and 
are  discussed  in  relation  to  their  subsequent  numerical 
evaluation.  By  determining  eigenvalues  of  equivalent 
depth  for  a  given  frequency  of  oscillation  instead  of  the 
other  way  around  the  problem  becomes  a  straightforward 
one  of  matrix  diagonalization..  If  solutions  are  based 
on  normalized  Legendre  polynomials  the  matrix  is  symmetric 
A  method  of  evaluating  the  related  wind  functions  is 
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1.  Introduction 

Laplace's  tidal  equation  governs  the  latitudinal 
variation  of  the  depth  of  a  thin  layer  of  uniform  fluid 
on  a  rotating  sphere  that  is  executing  small  periodic 
oscillations  in  longitude  <^>  and  time  t,  •  Non-zero 
solutions  are  found  for  discrete  values  (eigenvalues) 
of  the  frequency  of  oscillation  <r  :  alternatively,  if 
o'  is  given,  oscillations  are  found  for  eigenvalues  of 
the  undisturbed  depth  ft  . 

Laplace's  tidal  equation  also  arises  in  the  treatment 
of  oscillations  of  a  stratified  and  compressible  atmosphere 
h  being  replaced  by  a  constant  of  separation,  which  is 
usually  dinensionalized  to  a  certain  unit  of  length  and 
termed  an  equivalent  depth.  In  this  case  the  equation 
governs  the  latitudinal  variation  of  pressure  (and  also  of 
temperature  and  vertical  velocity)  at  any  given  height  and 
h  has  no  direct  physical  interpretation,  being  capable 
under  certain  conditions  of  assuming  negative  values. 

A  review  of  the  theory  of  atmospheric  tides  has  been 
presented  by  Chapman  and  Lindzen  (1970).  Various 
extensions  of  the  theory  have  since  been  made  to  damped 
oscillations  in  which  case  complex  values  of  CT  and  h 
are  introduced.  (Semenovskiy,  1971;  Volland,  197^a,b; 
Ishimine,  1977)* 

Work  on  the  solution  of  Laplace's  tidal  equation  was 


pioneered  by  Hough  (1897 »  I89Q)  in  terms  of  associated 
Legendre  polynomials  leading  to  a  third-order  recurrence 
relation  between  the  series  coefficients.  The  derivation 
of  this  relation  has  been  presented  on  several  occasions 
and  will  not  be  repeated  here*,  the  relation  is  to  be  found 
in  various  forms  which  have  led  to  different  numerical 
procedures  for  evaluating  the  eigenvalues.  This  paper 
x'eviews  the  various  theoretical  results  with  the  objective 
of  showing  their  relationship  and  providing  an  understanding 
relevant  to  their  numerical  analysis. 

2.  Laplace’s  tidal  equation 

Periodic  variations  in  longitude  and  time  have  been 
variously  formulated  in  terms  of  either  4,  or 

1  (^iji  —  ft ) 

4.  '  .  Westward  and  eastward  progressions  of  phase 

may  then  be  treated  by  assigning  both  positive  and  negative 
values  to  one  of  the  quantities  ,  s' .  We  limit  the 
discussion  to  two  forms  of  representation: 

(i)  -a.  ,  where  s  *  0,  +1,  +2,  ...  and  <r  is  a 

constant  positive  frequency,  and 
(ii;  jl  ,  where  s  is  a  non-negative  integer  and 

a  constant  non-zero  frequency. 


Form  (i)  was  adopted  by  Hough  (1897)  and  appears  more 
frequently  in  the  literature  than  form  (ii):  it  is  found 
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in  the  review  of  Chapman  and  Lindzen  (1970).  Form  (ii) 
has  been  adopted  by  Longuet-Higgins  (1968),  Volland  (1974a) 
and  Volland  and  hayr  (1977)*  With  (i)  a  westward 
(eastward)  progression  of  phase  has  s  positive  (negative); 
and  with  (ii)  a  westward  (eastward)  progression  of  phase 
has  <r  negative  (positive).  Any  results  obtained  in  one 
notation  may  be  readily  altered  to  apply  to  the  other 
according  to  the  scheme  in  Table  1. 

Table  1  s  Ranges  of  values  of  for 

westward  and  eastward  travelline:  waves. 

Direction  of  Form  (i)  Form  (ii) 

travel  of  wave 

Westward  r  >  o  ^>o  <r<o 

Eastward  ^<ot  <r  y  o  <s">o 


With  (i),  Laplace's  tidal  equation  may  be  written 


as 
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9  is  colatitude,  co  the  Earth’s  rotation  rate,  a  the 
Earth's  radius  and  g  the  acceleration  due  to  gravity. 
Table  1  and  (3)  show  that  with  form  (ii)  the  only  change 
required  to  (1)  is  the  replacement  of  s/f  by  -s/f  in  the 
centre  term  without  regard  to  the  direction  of  travel  of 
the  wave. 


The  recurrence  relation  with  form  (i)  periodic  terms 
For  given  s  (>0)  we  write 


0  =  1 


Substitution  of  (4)  into  (1)  leads  after  considerable 
reduction  to 


where 
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Equations  (5)  to  (8)  were  first  given  by  Hough  (1898). 

For  given  s  «0),  (1)  shows  that  ©  may  be  obtained 
by  replacing  s  by  -s  (>0)  and  f  by  -f  wherever  they  occur 
in  (6)  to  (8),  as  (1)  is  unaltered  by  these  changes  of 
sign.  For  s  =0,  M  is  infinite,  CQ  =  0  and  0  is 
based  on  • 

The  matrix  of  coefficients  of  (5)  has  the  form 

M  -  V  I  -  0  0 to) 

where  I  denotes  the  unit  matrix.  On  truncating  h  after 
the  n  th  row  and  column,  standard  computer  routines  are 
available  for  obtaining  n  eigenvalues  of  \  .  numerical 

accuracy  may  be  investigated  by  repeating  the  calculation 
with  different  values  of  n. 

In  Hough's  work,  fluid  depth  h  and  hence  X  were 
taken  as  given  and  eigenvalues  of  f  and  hence  <r  were 
sought.  Hough  replaced  (5>)  by  two  recurrence  relations 
involving  Cp  and  a  set  of  auxiliary  constants  Dr.  Then 


by  eliminating  Cr  and  Dr  an  equation  in  the  form  of  a 
continued  fraction  was  obtained  having  an  infinite  number 
of  roots  for  the  frequencies  of  equatorially  symmetric  modes. 
In  a  like  manner  another  continued  fraction  was  obtained 
for  asymmetric  nodes.  As  the  numerators  of  the  fractions 
were  independent  of  f,  only  the  denominators  needed  to  be 
re-determined  at  each  successive  iteration,  which  was 
particularly  advantageous  for  hand  calculations. 

If  is  given  and  eigenvalues  for  X  (and  hence  h) 
are  sought,  the  continued  fraction  iteration  can  still  be 
followed  but  appears  less  attractive  than  the  matrix 
formulation  (10)  as  starting  values  need  to  be  taken  for 
each  solution  whereas  diagonalization  of  M  yields  n 
eigenvalues  simultaneously.  In  recent  times,  eigenvalue 
solutions  by  the  continued  fraction  method  have  been 
evaluated  for  symmetric  modes  by  Kato  (1966)  and  Ishimine 
(1977). 

Chapman  and  Lindzen  (1970)  present  equation  (5)  and  give 
Hough's  two  sets  of  recurrence  relations  in  Cr  and  Dp  as 
the  basis  for  finding  eigenvalues  of  }\  and  hence  h. 
Procedurally ,  the  method  appears  less  attractive  than 
matrix  diagonalization  as  an  array  size  of  (2n  x  2n) 
instead  of  (n  x  n)  i3  needed  for  the  calculation  of  n 
eigenvalues  as  occurs  only  in  alternate  diagonal 
elements.  Also  as  the  form  of  (10)  is  lost,  an  iterative 
process  needs  to  be  followed  to  reach  each  root. 


4.  The  recurrence  relation  with  form  (ii)  periodic  terms 

Longuet-Higgins  (1968),  following  Love  (1913),  introduced 


functions  analogous  to  velocity  potential  and  stream  function 

$  =  S.  A-- p-  4'  o*) 


=  £  C '5r?rV> 


where  s  ^0.  The  recurrence  relations 
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were  obtained,  where 
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The  coefficient  of  Ar_^  in  (13)  should  read  -qr_^  in 
Longuet-Higgins  (1968)  instead  of  _9r+^»  and  in  equation 
(3«19)  of  that  paper  the  sign  of  the  term  in  B®  needs  to 
be  changed.  The  same  equation  has  been  reproduced  in 
equation  (28)  of  Volland  and  liayr  (197?)  without  correction 
and  the  same  error  of  sign  appears  in  equation  (2.12)  of 


Volland  (1974-a)  where  the  detei-minant  of  the  coefficients 
in  (13) ,  namely 


K*  |V*+I  0  0  *  *  * 

L/4+t  0  ,  .  . 

0  ^3  ‘  '  ' 


(Iff) 


appears  without  the  two  negative  signs.  On  eliminating 
from  (13)  we  obtain 


In  the  case  of  a  uniform  fluid,  the  dependent  variable 
in  (1)  may  be  taken  as  the  fluid  depth  and  hence  as  the 
rate  of  change  of  fluid  depth  with  time.  By  mass  continuity 
the  rate  of  decrease  of  fluid  depth  equals  the  divergence 
of  fluid  from  a  unit  column,  and  this  is  proportional  to 

,  where  SJl  denotes  the  horizontal  Laplacian  operator 
on  a  unit  sphere.  Hence  0  in  (1)  (with  s/f  replaced  by 


m 


-s/f)  is  such  that 


In  the  case  of  a  compressible  atmosphere,  the  proportionality 
is  height  dependent.  From  (4),  (11)  and  the  relation 

v  0«) 

r  • 

it  follows  that 

Ar  cr/rc^o  c^> 

On  substituting  (14)  and  (19)  into  (16)  a  relation  in 
is  obtained  which  agrees  with  (5)  on  changing  s/f  to  -s/f 
to  accord  with  the  use  of  form  (ii). 

The  matrix  of  coefficients  in  (13)  is  seen  by  (14)  to 
separate  into  terms  in  f  and  \/f  having  the  form 
(Longuet-Higgins,  1968) 

tv/f)  j  +  c  -  /  i  <?*) 

where  J  and  C  are  readily  constructed  from  (14).  For 
given  Vf  =  cQ,  say,  eigenvalues  of  f  may  be  obtained  by 
matrix  diagonalization  of  cQJ  +  C.  This  procedure  does 
not  directly  solve  the  problem  of  finding  eigenvalues  of 
f  for  given  A  ,  but  provides  a  sequence  of  pairs  of 
values  (fi,  Ai)  such  that  Aj/f^  -  c0*  By  taking  a 
range  of  values  of  cQ,  Longuet-Higgins  (1968)  was  able 
to  plot  curves  of  f  =  f(A),  each  curve  relating  to  a 
particular  mode.  Alternatively  if  a  range  of  values 
of  f  were  taken,  the  same  set  of  curves  in  the  form 
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X  »  X(f)  could  be  obtained  by  diagonalizing  M  of  equation 
(10)  (with  s/f  changed  to  -s/f).  The  latter  procedure 
was  followed  in  principle  by  Volland  (1974a ,b)  in  an 
investigation  of  the  effect  on  X  (and  hence  on  h)  of 
varying  the  imaginary  part  of  a  complex  f  for  particular 
modes.  A  matrix  of  coefficients  in  the  form  (20)  was 
adopted  and  on  noting  that  only  alternate  members  of  the 
diagonal  elements  Kg,  Kg+2»«*»  contained  X,  the 

determinant  was  condensed  to  one  of  half  the  order  by 
operations  between  columns.  The  effect  of  this  reduction 
was  the  same  as  eliminating  from  (13)  and  using  (16). 

5.  The  use  of  normalized  P„  „ 

_ r,s 

The  results  reviewed  above  have  been  based  on  the 
associated  Legendre  polynomials 

p^  -  0-h')  —r  (?•«) 

r  i  ..a 


which  are  related  to  the  normalized  form  P„  „ 

r,s 


by 
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Writing 
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and  comparing  with  (4)  after  substituting  from  (22)  we 


have 


(it 

2.  (r+^  }  ■ 


(ft*.) 


Substituting  (24)  into  (5)  and  multiplying  through 

L 

by  ^2(r  +  s)!/(2r  +  l)(r  -  s)^]x,  we  obtain 

Lr-aar-*  -v  CM,-V)ar  *  |___  o.^  =  O  taO 

where 

i 

K  (Xrn  3)  [t.ar-*  i)  (5r  4-y^J  v  -  (r-n)(r* 

Equations  (25)  and  (26)  have  been  obtained  by  Dikii  (1965) 
although  the  factor  (2r  +  5)  in.  the  denominator  of  (26) 
was  misprinted  as  (2r  +  s).  A  typographical  error  also 
appears  in  the  recurrence  relations  which  precede  the 
derivation  of  (25) *  where  the  factor  (s/f  -  n  +  l) 
appearing  in  the  coefficient  of  bR_^  should  read 
(s/f  -  n  -  1).  Dikii  (1965)  obtained  the  expression 
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Equation  (27)  may  be  shown  to  be  the  same  relation  as 
(7)  by  separating  it  into  its  component  partial  fractions 
with  denominators  (s/f  *  r),  (s/f  -  r  -  1),  [s/f  -  r(r  -  l)~] 
and  £s/f  -  (r  +  l)(r  +  2)]  and  then  combining  terms  in 
(s/f  +  r)  and  (s/f  -  r  -  1). 

From  (25)  the  matrix  of  coefficients  is  obtained  as 
F  -  X  I ,  where 


F 


0  0  0 

0  %  0  L*,,  0  ... 

L"  0  0  .  .  . 

0  0  0  .  .  * 

0  0  l**  o  .  •  • 

'V4'4- 

*  •  •  • 


(a») 


As  F  is  symmetric,  it  may  be  shown  (Jones,  1970)  that 

the  eigenvectors  are  orthogonal  and  that,  on  account  of 

the  orthogonality  of  P_  the  eigenfunctions  (Hough 

r,s 

functions)  are  orthogonal. 

Normalized  spherical  harmonics  were  adopted  by  Jones 
(1970)  in  the  development  of  a  general  theory  of  atmospheric 
oscillations  by  expanding  in  terms  of  vector  harmonics  to 
obtain  a  matrix  formulation  of  the  equations  of  motion. 

On  approximating  to  the  case  of  classical  tidal  theory, 
mode  decoupling  is  achieved  by  transforming  to  basic  field 
variables  that  diagonalize  a  certain  matrix  of  infinite 
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order  whose  inverse  is 


L’  [A  -  fTA-TlL" 


where  T  has  non-zero  elements 
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L  and  A  are  diagonal  matrices  whose  diagonal  elements 
have  r  *  8,  s  +  1,...  • 

On  multiplying  out  (29) ,  it"  is  found  that 

H"'  -  p  F  (■«) 

—1  p 

The  eigenvalues  of  H  are  therefore  those  of  F  times  f  , 
i.e*  f2/ )vi  or  a2cr2/ehi  by  (3)  and  (9). 

For  numerical  work,  the  order  of  arrays  may  be  halved 
by  treating  (23)  as  two  sets  of  equations  whose  coefficients 
form  matrices  FQ  -  "X  I  and  F^  -  Xl,  where 

t\.  .  i  w 
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The  first  set  leads  to  eigenvectors  for  (ag,  as+2»»**  ) 
and  hence  to  eigenfunctions  by  (25)  that  are  equatorially 

symmetric  being  based  on  P3,S-  Fs*2 .  The  second  set 

leads  to  eigenvectors  for  (ag+1,  as+2’***  ^  which  are  based 
on  Pg+^  a»  ^s+3  g » •  •  •  and  hence  to  the  asymmetric  solutions 
(Dikii,  1965). 


6.  ^valuation  of  eigenfunctions  and  wind  functions 
before  calculating  (h)  from  (25),  the  vector 
(as,  0,  ag+2,  0,...  )  is  normalized  by  dividing  it  by 

(ag  +  as+2  +  ***  ^  *  an^  ^as+l’  as+5*  0,.».  )  i® 
treated  likewise.  ©  is  then  determined  apart  from  its 

sign,  which  is  arbitrary.  Signs  may  be  chosen  so  that 
symmetric  functions  are  positive  at  the  equator  and 
asymmetric  functions  arc  increasing  with  latitude  at  the 
equator. 

The  horizontal  wind  components  of  a  tidal  oscillation 
depend  on  horizontal  gradients  of  the  pressure  field  and 
latitudinally  depend  on  functions  derived  from  ©  •  With 
form  (i),  these  functions  are  (Chapman  and  Lindzen,  1970) 
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for  the  eastward  and  northward  components  respectively. 

Series  expansions  of  0O  ,  0^  in  terms  of  Pr  g  are 
not  obtainable,  but  developments  are  possible  for 


Su  -  6  0O 

Prom  (36)  and  (37)  we  have 


5.,  -  -^<9  0, 


(?*) 


f  5",  -  fS*  -  (W|)0  -  o 
hS„  -  fSv  -T0  '  • 


where  3)  -  (i-  •  For  s  ^  0,  we  write 


$U  B  ^  vr  . 


and  note  that  (39)  is  unaltered  for  8  >  0  if  written  as 


f  s„  -  p5v  -  t'-'/f)© 


-  O 


L^) 


Then  by  (23)  ♦  (41)  and  the  relations 
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where 

ir  ^  [Crt^yCA.r'-.)]i  t-k*) 

the  coefficients  of  P_  ,  in  (40)  and  (42)  give 

y  s 

*  i  vr-t  “  ^r-i  'w-r-a  ~  Cr"4^r-i  ar-i+ (r+')dror  (Ub) 

&rVr  ~  j  *r„{  -  ^r.(  trr.2  -  CM,/p  ar-,  (A7^ 


By  successively  putting  r  =  s,  s  +  1,...  and  taking 
ar  »  up  »  vr  *  0  for  r  <  sf  (46)and  (47)  enable 

(us»  °»  us+2’  °»***  )  and  (°»  vs+i»  °»  vs+5» • • •  )  to  be 
calculated  for  symmetric  Hough  functions  for  which 

as+l  =  as+3  *  •••  "  °»  and  likewise  (0,  ug+1,  0,  ug+5,...  ) 
and  (v  t  0,  vg+2»  0,...  )  for  asymmetric  Hough  functions 
for  which  ag  ■  ag+2  ~  “  0*  0y  »  can  then  be 

obtained  from  1,33;  and  (41)  except  for  the  end  points 
B  *  0,  ic  . 

To  calculate  when  s  <  0,  it  is  only 

necessary  to  replace  s  by  -s  (>0),  and  f  by  -f,  as  it 
follows  from  (40)  and  (42)  that  the  above  procedure  leads 


to  the  evaluation  of  0^  and  -  0^.  If  f  <  1,  (36)  and 
(37)  become  indeterminate  at  |*.»  f,  whereas  no  such 
difficulty  arises  with  the  above  procedure. 

7«  Discussion 

The  recurrence  relation  (5),  which  is  that  originally 
derived  by  Hough,  has  a  matrix  of  coefficients  of  the 
form  of  equation  (10)  and  standard  computing  routines 
are  available  for  the  determination  of  its  eigenvalues 
and  eigenvectors.  If  eigenvalues  of  equivalent  depths  are 
required  for  a  given  frequency  of  oscillation  the  results 
are  obtainable  in  one  step,  but  if  eigenvalues  of  frequency 
are  required  for  a  given  depth  an  approximating  iteration 
is  necessary.  The  latter  case  arises  in  a  study  of  the 
free  oscillations  of  the  atmosphere  as  equivalent  depth 
is  then  determined  as  an  eigenvalue  of  the  vertical 
structure:  it  was  ilso  the  problem  that  concerned  Hough 
and  in  devising  an  iterative  procedure  he  undertook 
further  analytical  developments  of  the  recurrence  relation. 
These  developments  are  not  however  essential  to  the 
eigenvalue  analysis  and  it  it  doubtful  whether  they  offer 
any  advantage  over  the  matrix  formulation  with  present-day 
computing  facilities. 

The  recurrence  relation  recommended  for  use  is 
equation  (23)  which  was  given  by  Dikii  and  corresponds 


-  20  - 


to  solutions  based  on  normalized  associated  Legendre 
polynomials.  Apart  from  th£  better  numerical  conditioning 
that  is  to  be  expected  with  normalization,  there  is  the 
simplification  of  a  symmetric  matrix. 
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